Effects of density and parametrization on scattering observables 
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We calculate the density distribution of protons and neutrons for 40>'*2,44,48^^ frame- work of relativistic 

mean field (RMF) theory with NL3 and G2 parameter sets. The microscopic proton-nucleus optical potential 
for p -I-*" Ca system is evaluted from Dirac NN-scattering amplitude and the density of the target nucleus using 
Relativistic-Love-Franey and McNeil- Ray- Wallace parametrizations. Then we estimate the scattering observables, 
such as elastic differential scattering cross-section, analysing power and the spin observables with relativistic 
impulse approximation. We compare the results with the experimental data for some selective cases and found 
that the use of density as well as the scattering matrix parametrization is crucial for the theoretical prediction. 



Explaining the nuclear structure by taking the 
tool of nuclear reaction is one of the most curious 
and challenging solution for Nuclear Physics both 
in theory and laboratory. So far the elastic scat- 
tering reaction of Neucleon- Nucleus is more inter- 
esting than that of Nucleus-Nucleus at laboratory 
energy Eiab 1000 MeV. The Neucleon-Nucleus 
interaction provides a fruitful source to determine 
the nuclear structure and a clear path toward the 
formation of exotic nuclei in laboratory. One of 
the theoritical method to study such type of re- 
action is the Relativistic Impulse Approximation 
(RIA). In a wide range of energy interval, the 
conventional impulse approximation |1|2) repro- 
duces quantitatively the main features of quasi- 
elastic scattering for medium mass nuclei [314] . 
The observables of the elastic scattering reaction 
not only depend on the energy of the incident 
particle but also on the kinematic parameter as 
well as the density discributions of the target nu- 
cleus. In the present letter, our motivation is to 
calculate the nucleon-nucleus elastic differential 
scattering cross-section (^) and other quanti- 
ties, like optical potential (Uopt), analysing power 
(Ay) and spin observables (Q— value) taking in- 
put as relativistic mean field (RMF) and recently 
proposed effective field theory motivated rela- 
tivistic mean field (E-RMF) density. The RMF 
and E-RMF densities are obtained from the most 



successful NL3 [5] and advanced G2 [6] param- 
eter sets, respectively. As representative cases, 
we used these target densities folded with the 
NN-aplitude of 1000 MeV energetic proton pro- 
jectile with Relativistic-Love-Franey (RLE) and 
McNeil-Ray- Wallace (MRW) parametrizations [7] 
for 40,42,44,48(^Q -^^ q^j. calculations. 

The RMF and E-RMF theories are well docu- 
mented |6I8I9] and for completeness we outline 
here very briefly the formalisms for finite nu- 
clei. The energy density functional of the E-RMF 
model for finite nuclei is written as |10|llj . 
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NN-scattering operater F can be written in terms 
of five complex functions (the five terms involves 
in the proton-proton pp and neutron-neutron pn 
AAR, (1) scattering) as follows: 



where the index a runs over all occupied states 
(/^^(r) of the positive energy spectrum, $ = 
5s'/'o(r), W EE g„Vb(r), R = gpbo{r) and A = 
eAo(r). The terms with A, /3s and take care 
of the effects related with the electromagnetic 
structure of the pion and the nucleon (see Ref. 

The energy density contains tensor cou- 
plings, and scalar- vector and vector-vector meson 
interactions, in addition to the standard scalar 
self interactions K3 and K4. Thus, the E-RMF for- 
malism can be interpreted as a covariant formu- 
lation of density functional theory as it contains 
all the higher order terms in the Lagrangian, ob- 
tained by expanding it in powers of the meson 
fields. The terms in the Lagrangian are kept fi- 
nite by adjusting the parameters. Further insight 
into the concepts of the E-RMF model can be 
obtained from Ref. [11 . It may be noted that 
the standard RMF Lagrangian is obtained from 
that of the E-RMF by ignoring the vector- vector 
and scalar-vector cross interactions, and hence 
does not need a separate discussion. In each of 
the two formalisms (E-RMF and RMF), the set 
of coupled equations are solved numerically by a 
self-consistent iteration method and the baryon, 
scalar, isovector, proton, neutron and tensor den- 
sities are calculated. 

The numerical procedure of calculation and the 
detailed equations for the ground state properties 
of finite nuclei, we refere the reader to Refs. |9|8j . 
The densities obtained from RMF (NL3) [5] and 
E-RMF (G2) are used for folding with the NN- 
sacttering amplitude at Eiah = lOOOAfeF, which 
gives the proton-nucleus complex optical poten- 
tial for RMF and E-RMF formalisms. RIA in- 
volves mainly two steps |12|13j of calculations for 
the evaluation of the NN-scattering amplitude. 
In this case, five Lorentz covariant function f7 
multiply with the so called Fermi invariant Dirac 
matrix (NN-scattering amplitudes). This NN- 
amplitudes are folded with the target densities 
of protons and neutrons to produced a first or- 
der complex optical potential f/opt ■ The invariant 



^ 7(0)7(1) 



+-^^<o')^(i)M'v +-^^7(^0)7(%7(\)7(i)M. (2) 

where (0) and (1) are the incident and struck nu- 
cleons respectively. The amplitude for each 
is a complex function of the Lorentz invariants T 
and S with E = Eiab and q is the four momentum. 
We recommend the redears for detail expressions 
to Refs. |14|15|16|17|18I19|20|21|22] . Then the 
Dirac optical potential Uopt{q,E) can be written 
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where F is the scattering operator, p is the mo- 
mentum of the projectile in the nucleon-nucleus 
center of mass frame, |^) is the nuclear ground 
state wave function for A-particle. Finally using 
the Numerov algorithm the obtained wave func- 
tion is matched with the coulomb scattering solu- 
tion for a boundary condition at r — > cxd and we 
get the scattering observables from the scattering 
amplitude, which are defined as: 
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Now we present our calculated results of neu- 
trons and protons density distribution obtained 
from the RMF and E-RMF formalisms [8]. Then 
we evaluate the scattering observables using these 
densities in the relativistic impulse approxima- 
tion, which involves the following two steps: 
in the first step we generate the complex NN- 
interaction from the Lorentz invariant matrix 
F^{q,E) as defined in Eq. (2). Then the in- 
teraction is folded with the ground state target 
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nuclear density for both the RLF and MRW pa- 
rameters [7; separately and obtained the nucleon- 
nucleus complex optical potential UoptiQ, E) for 
the parametrisations. It is to be noted that pair- 
ing interaction is taken care using the Pauli block- 
ing approximation. In the second step, we solve 
the wave fonction of the scattering state utilising 
the optical potential prepared in the first step by 
well known Numerov algorithm |23| . The result 
approxumated with the non-relativistic Coulomb 
scattering for a longer range of radial component 
which results the scattering amplitude and other 
observables [211 ■ thr present paper we calcu- 
late the density distribution of protons and neu- 
trons for 40,42,44,48^^^ jy^^g Q2 parameter 

sets. From the density we evalute the optical po- 
tential and other scattering observables and some 
representative cases are presented in Figures 1—3. 
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Figure 1. (upper panel): The neutrons and 
protons density distribution for ^^Ca with NL3 
and G2 parameter sets, (lower panel) (a) the 
Dirac optical potential for p -f**^ C'a system us- 
ing RMF (NL3) and E-RMF (G2) densities with 
RLF parametrisation, (b) same as (a), but for 
MRW parametrisation. The projectile proton with 
Eiab = 1000 MeV IS taken. 



In Fig. 1, the protons and neutrons density 
distribution for ■*"Ca using NL3 and G2 parame- 
ter sets (upper panel) and the optical potential 
obtained with RLF and MRW parametrisation 
for p -)-^° Ca at 1000 MeV proton energy (lower 
panel) are shown. From the figure, it is noticed 
that, there is no significant difference in desities 
for RMF and E-RMF parameter sets. However, 
a careful inspection shows a small enhancement 
in central density (0-1.6 fm) for NL3 set. On the 
otherhand the densities obtained from G2 elon- 
gated to a larger distance towards the tail part 
of the density distribution. As the optical poten- 
tial is a complex function which constitute both 
real and imaginary part for both scalar and vec- 
tor, we have displyed those values in the lower 
panel of Fig. 1. Unlike to the (upper panel) of 
protons and neutrons density distribution, here 
we find a large difference of Uopt{q,E) between 
the RLF and MRW parametrisation. Further, the 
Uopt{q, E) value of either RLF or MRW differs sig- 
nificantly depending on the NL3 or G2 force pa- 
rameters. That means, the optical potential not 
only sensitive to RLF or MRW but also to the use 
of NL3 or G2 parameter sets. Investigating the 
figure it is clear that, the extrimum magnitude of 
real and imaginary part of the scalar potential are 
-442.2 and 113.6 MeV for RLF (G2) and -372.4 
and 109.1 MeV for RLF (NL3). The same values 
for the MRW parametrisation are -219.8 and 32.8 
MeV with G2 and -175.1 and 33.2 MeV with NL3 
sets. In case of vector potential, the extremum 
values for real and imaginary parts are 361.3 and 
-179.2 MeV for RLF (G2) and 279.2 and -164.8 
MeV for RLF (NL3) but with MRW parametrisa- 
tion these are appeared at 128.1 and -87.4 MeV in 
G2 and 99.2 and -76.6 MeV in NL3. From these 
large variation in magnitude of scalar and vector 
potentials, it is clear that the predicted results not 
only depend on the input target density, but also 
highly sensitive with the kinematic of the reaction 
dynamics. A further analysis of the results for 
the optical potential with NL3 and G2, it suggest 
that the Uopt value extends for a larger distance in 
NL3 than G2. For example, with RLF the central 
part of Uopt with G2 is more expanded than with 
NL3 and ended at r ^ 6/m, whereas the optical 
potential persists till r ^ 8/m in NL3. Similar 
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situation is also valid in MRW parametrisation. 
This nature of the potential suggests the appli- 
cability of NL3 over G2 force parameter. This is 
because in case of NL3 the soft-core interaction 
between the projectile and the target nucleon is 
more effective. 
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Figure 2. The elastic differential scattering cross- 
section (^) as a function of scattering angle 

Ocm(deg) for 40,42,44,48(-(j ygj^^ ^0^^ Jl^p 

MRW parametrisations. The value of ^ is 
shown for RMF (NL3) and E-RMF (G2) densi- 
ties. 



In Fig. 2., we have plotted the elastic scatter- 
ing cross-section of the proton with 40,42,44,48^^^ 
at laboratory energy Eiab =1000 MeV using both 
densities obtained in the NL3 and G2 parameter 
sets with RLF and MRW parametrisations. The 
experimental data |25j are also given for compar- 
ison. It is reported in Refs. |7|26j the superiority 
of RLF over MRW for lower energy [Eiab < 400 
MeV), however the MRW shows better results at 
energy Eiab > 400 MeV. In the present case, our 
incident energy is 1000 MeV which matches bet- 
ter (MRW) with experimental values. This is con- 
sistent with the optical potential also (see Fig. 
1). From the differential cross-section for both 
NL3 and G2 densities with MRW parametriza- 



tion, it is clearly seen that ^ with NL3 desity is 
more closer to experimental data which insist not 
only the importance of parametrization (RLF or 
MRW) but also to choose proper density input 
for the reaction dynamics. Analysing the elastic 
differential cross-section along the isotopic chain 
of Ca from A=40 to 48, the calculated results im- 
prove with increasing mass number of the target. 
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Figure 3. (a) The calculated values of analysing 
power Ay as a function of scattering an- 
gle 6cm (deg) for "^^Ca (b) The spin observ- 
able Q— value as a function of scattering angle 
Qcm(deg) for *^Ca. In both (a) and (b), the RLF 
and MRW parametrisations are used with RMF 
(NL3) and E-RMF (G2) densities. 



The analysing power for p -l-^*' Ca composite 
system is calculated from the general formulae 
given in eqns. (4) and (5) and are shown in Fig. 
3 with RLF and MRW. The Ay and Q-values ob- 
tained by NL3 and G2 sets almost matches with 
each other both in RLF and MRW. But if we com- 
pare the results with RLF and MRW it differs sig- 
nificantly. Again, we get a small oscillation of Ay 
and Q in G2 set with increasing scattering angle 
0^,-^ for RLF which does not appear in NL3 set. 
There is a rotation of Q— value from positive to 
negative direction when we calculate with MRW 
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parametrization, which does not appear in case 
of RLF parametrization. This rotation shows a 
shining path towards the formation of exotic nu- 
clei in the laboratory. 

In summary, we calculate the density distri- 
bution of protons and neutrons for 4o,42,44,48(-^^ 
by using RMF (NL3) and E-RMF (G2) param- 
eter sets. We found similar density distribution 
for protons and neutrons in both the sets with a 
small difference at the central region. This small 
difference in densities make a significant influence 
in the prediction of optical potential, elastic dif- 
ferential cross-section, analysing power and the 
spin observable for p + Ca systems. The effect of 
kinematic parameters for reaction dynamics, RLF 
and MRW, are also highly sensitive to the pre- 
dicted results. That means, the differential scat- 
tering cross-section and scattering observables are 
highly depent on the input density and the choice 
of parametrisation. 
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